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ABSTRACT: 

Accurate quantum Monte Carlo calculations of ground and low-lying excited states of light p-shell 
nuclei are now possible for realistic nuclear Hamiltonians that fit nucleon-nucleon scattering data. 
At present, results for more than 30 different (J n ;T) states, plus isobaric analogs, in A < 8 nuclei 
have been obtained with an excellent reproduction of the experimental energy spectrum. These 
microscopic calculations show that nuclear structure, including both single-particle and clustering 
aspects, can be explained starting from elementary two- and three-nucleon interactions. Various 
density and momentum distributions, electromagnetic form factors, and spectroscopic factors have 
also been computed, as well as electroweak capture reactions of astrophysical interest. 
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1 INTRODUCTION 



A major goal of nuclear physics is to understand the stability, structure, and re- 
actions of nuclei as a consequence of the interactions between individual nucleons. 
This is an extremely challenging many-body problem, exacerbated by the fact that 
we do not know in detail what the interactions are. Quantum chromodynamics, 
the fundamental theory for strong interactions, is so difficult to solve in the non- 
perturbative regime of low-energy nuclear physics that no one is currently able 
to quantitatively describe the interactions between two nucleons from first princi- 
ples. However, we do have a tremendous amount of experimental data on nucleon- 
nucleon (WW) scattering and we can construct accurate representations of AW 
interactions by means of two-body potentials. These potentials are complicated, 
depending on the relative positions, spins, isospins, and momenta of the nucle- 
ons, which makes the problem of finding accurate quantum-mechanical solutions 
of many-nucleon bound and scattering states a demanding task. Nevertheless, a 
variety of precise calculations have been carried out for the three- and four-nucleon 
systems which reveal an additional complication: two-nucleon forces alone are in- 
adequate to quantitatively explain either the bound-state properties or a variety 
of Nd and Nt scattering data. This is not surprising, because meson-exchange 
theory and the composite nature of nucleons indicate that many-nucleon forces, 
at least three-nucleon forces, are important. However, these are even more diffi- 
cult to construct from first principles. In like manner, meson-exchange introduces 
two-body charge and current operators, and calculations of electroweak transitions 
show that these are also necessary to achieve agreement with data. 

Despite these difficulties, tremendous progress has been made in the past decade, 
both in the characterization of nuclear forces and currents, and in the development 
of accurate many-body techniques for evaluating them. A multienergy partial- 
wave analyses of elastic AW scattering data below = 350 MeV was produced 
by the Nijmegen group in 1993 (|l]), which demonstrated that over 4300 data points 
could be fit with a x 2 P er datum ~ 1. This was achieved by careful attention to 
the known long-range part of the AW interaction and a rigorous winnowing of 
inconsistent data. Their work inspired the construction of a number of new WW 
potential models which gave comparably good fits to the database. 

Concurrently, several essentially exact methods have been developed for the 
study of few-nucleon (A = 3,4) systems with such realistic WW potentials. The 
first accurate, converged calculations of three-nucleon (WWW) bound states were 
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obtained by configuration-space Faddeev methods in 1985 (Q), while the first accu- 
rate 4 He bound states were obtained by the Green's function Monte Carlo method 
in 1987 (|3|). Other methods of comparable accuracy have been developed since then 
for both A = 3,4 bound states and for NNN scattering. A comprehensive review 
of many aspects of the few-nucleon problem can be found in the article by Carlson 
and Schiavilla (||). The result of all the progress in the last decade is that we 
have obtained an excellent understanding of few-nucleon systems, including both 
structure and reactions, although a few unanswered problems remain. 

Until recently, nuclei larger than A = 4 have been described within the frame- 
work of the nuclear shell model or mean- field pictures such as Skyrme models. 
The classic description of p-shell (5 < A < 16) nuclei was developed by Cohen 
and Kurath in the 1960s In that work, a large number of experimental levels 
were fit in terms of a relatively small number of two-body matrix elements, but no 
direct connection was made to the underlying NN force. Beginning with the work 
of Kuo and Brown (|6|), approximate methods for relating the shell- model matrix 
elements to the underlying forces were developed and applied to even larger nuclei. 
However, in most calculations, only the valence nucleons, i.e., the nucleons in the 
last unfilled shell, are treated as active, while the core is inert. Recently the first 
"no-core" shell model calculations, in which all the nucleons are active, have been 
made in the p-shell with effective interactions that have been derived by a G-matrix 
procedure from a realistic underlying AW force (|^). These calculations agree with 
the other exact methods for the ^4 = 3,4 nuclei, but fully-converged calculations 
for larger systems have not yet been obtained. 

The present article focuses on recent developments in quantum Monte Carlo 
methods that are making the light p-shell nuclei accessible at a level of accuracy 
close to that obtained for the s-shell (A < 4) nuclei. The quantum Monte Carlo 
methods include variational Monte Carlo (VMC) and Green's function Monte Carlo 
(GFMC) methods. The VMC is an approximate method (||,^,[l(],[ll]) that is used as 
a starting point for the more accurate GFMC calculations. Because of its relative 
simplicity, the VMC has also been used to make first studies of various reactions, 
such as the nuclear response to electron (|T^,|l3|) or pion scattering (ll4|), and for 
electroweak capture reactions of astrophysical interest ( |i~5| , |i~6| ) . The GFMC is exact 
in the sense that binding energies can be calculated with an accuracy of better 
than 2%; it has been used to compute ground states and low-lying excitations 
of all the A < 8 nuclei ^3, 17, 18,|l9). These studies of light p-shell nuclei have 
shown that nuclear structure can be predicted from the elementary nuclear forces. 
As presently implemented, both the VMC and GFMC methods require significant 
additional computer resources for every nucleon added, so they probably will be 
limited to systems of A < 12 nuclei. A new method, whose accuracy lies between 
that of the VMC and GFMC methods, is auxiliary-field diffusion Monte Carlo 
(AFDMC); it has been implemented for large pure neutron systems (20) and with 
further development should be capable of handling nuclei above A = 12. 

Most of the calculations discussed here have been made using Hamiltonians con- 
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taining the Argonne Vig NN potential ([2~I|) alone or with one of the Urbana (|j|) or 
Illinois ( |2^ ) series of NNN potentials. Argonne Vig (AV18) is representative of the 
modern NN potentials that give accurate fits to scattering data, while the Urbana 
and Illinois models are modern NNN potentials based on meson-exchange inter- 
actions and fit to the binding energies of light nuclei. These models are described 
in some detail in Section 2 to show the complexity of modern nuclear forces, and 
the consequent challenge for many-body theory. The VMC method is presented in 
Section 3; it starts with the construction of a variational wave function of specified 
angular momentum, parity and isospin, ^'y(J 7r ;T), using products of two- and 
three-body correlation operators acting on a fully antisymmetrized set of one-body 
basis states coupled to specific quantum numbers. Metropolis Monte Carlo inte- 
gration ([23]) is used to evaluate (\Ev|.ff|\E r v'), giving an upper bound to the energy 
of the state. The GFMC method, described in Section 4, is a stochastic method 
that systematically improves on ^>y by projecting out excited state contamination 
using the Euclidean propagation *(r) = exp[— {H — Eq)t]^v- In the limit r — > 00, 
this leads to the exact 

Energy results for the ground states and many low-lying excited states of the 
A < 8 nuclei are presented in Section 5. With the new Hamiltonians, the excitation 
spectra are reproduced very well, including the relative stability of different nuclei 
and the splittings between excited states. Energy differences between isobaric 
analog states and isospin-mixing matrix elements have also been computed to study 
the charge- independence-breaking of the strong and electromagnetic forces. A 
variety of density and momentum distributions, electromagnetic moments and form 
factors, transition densities, spectroscopic factors, and intrinsic deformation are 
presented in Section 6. One of the chief benefits of this approach is that absolute 
normalizations can be computed, and effective charges so commonly needed in 
traditional shell-model calculations of electroweak transitions are not required. 

While the quantum Monte Carlo methods have been used primarily for station- 
ary states, it is also possible to study scattering states, as discussed in Section 7. 
Some first calculations of astrophysical radiative capture reactions in the p-shell 
are presented in Section 8. Both of these applications are in their early stages, but 
promise to be major fields of endeavor for future work. In Section 9 we discuss 
neutron drops, systems of interacting neutrons confined in an artificial external 
well, that can serve as reference points for the development of Skyrme interactions 
for use in very large neutron-rich nuclei. Finally, in Section 10, we consider the 
outlook for future developments. 

2 HAMILTONIAN AND CURRENTS 

The Hamiltonian used in this review use includes nonrelativistic kinetic energy, the 
AV18 NN potential (§l[) and either the Urbana IX @ or one of the new Illinois 
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i i<j i<j<k 

The kinetic energy operator has charge-independent (CI) and charge-symmetry- 
breaking (CSB) components, the latter due to the difference between proton and 
neutron masses, 

IS = K? 1 + K? 3B = (- + — )V? - " — , (2) 

where Tj is the isospin of nucleon i. 

The AV18 potential is written as a sum of electromagnetic and one-pion-exchange 
terms and a shorter-range phenomenological part, 

Vij = v 7. + V 7. + V R . ( 3 ) 

The electromagnetic terms include one- and two-photon-exchange Coulomb inter- 
action, vacuum polarization, Darwin-Foldy, and magnetic moment terms, all with 
appropriate proton and neutron form factors. The one-pion-exchange part of the 
potential includes the charge-dependent (CD) terms due to the difference in neutral 
and charged pion masses. It can be written as 

r2 I "' \ : i' 



wy = / l— i mc X ijTi ■ tj + , (4) 

where TJj = 3r z jT z j — Tj • tj is the isotensor operator and 

Xi i = 3 (^ij + ' ( 5 ) 

X l3 = I (XP- - X±) , (6) 
X™ = [Yimr^Vi ■ <x, + T(mr^)^] . (7) 

Here Y(x) and T(x) are the normal Yukawa and tensor functions: 

Y(z) = — (8) 
x 

T(x) = ( -I + - + l) Y (x) fr(r) . (9) 

with short-range cutoff functions £y(r) and ^T( r ) ; and the X ±,D are calculated 
with m = and m n o. The coupling constant in Eq.(||) is f 2 = 0.075 and the 
scaling mass m s = m^i . 

The remaining terms are of intermediate (two-pion-exchange) and short range, 
with some 40 adjustable parameters. The one-pion-exchange and the remaining 
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Figure 1: Three-body force Feynman diagrams: a is the Fujita-Miyazawa or two- 
pion P-wave, b is two-pion S-wave, c and d are three-pion rings with one A in 
intermediate states. 



phenomenological part of the potential can be written as a sum of eighteen opera- 
tors, 

v ij+ v S= E «p(^)ofi- (io) 

p=l,18 

The first fourteen are CI operators, 



r)P=i,i4 

ij 



[l,<Tj • 0-j,Sij,li S,L ,L (<Ti ■ <Tj), (L • S) ] <g) [l,Ti ■ Tj] , (11 

while the last four, 

[1, cTi ■ (Tj,Sij] <g> Tij , and (r zi + r zj ) , 



=15,18 



(12) 



are CD and CSB terms. The potential was fit directly to the Nijmegen NN scat- 
tering data base ([!]) containing 1787 pp and 2514 np data in the range 0-350 MeV, 
to the nn scattering length, and to the deuteron binding energy, with a x 2 P er 
datum of 1.09. 

While the CD and CSB terms are small, there is a clear need for their presence. 
The Nijmegen group studied many NN potentials from the 1980s and before, and 
found that potentials fit to np data in T = 1 states did not fit pp data well even 
after allowing for standard Coulomb effects, and vice versa (^). All five modern 
NN potentials, Reid93, Nijm I and II (H), CD Bonn (||), and AV18, which fit pp 
and np scattering data with a xVpoint near 1, have CD components. The CSB 
term is required to fit the difference between the pp and nn scattering lengths, and 
is consistent with the mass difference between 3 H and 3 He, as shown below. The 
identification of the proper CI, CD, and CSB NN force components is important in 
setting the proper baseline for the introduction of NNN forces, which are required 
to fit the 3 H and 3 He binding energies as a primary constraint. 

The Urbana series of three-nucleon potentials is written as a sum of two-pion- 
exchange P-wave and remaining shorter-range phenomenological terms, 



V, 



ijk 



(13) 
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The structure of the two-pion P-wave exchange term with an intermediate A ex- 
citation (Fig. |l]a) was originally written down by Fujita and Miyazawa (^); it can 
be expressed simply as 

= E (^{^.^H^ • Tj'Ti ■ Tk} + C&[X%, Xj k ] [ Ti • Tj , Tj ■ Tfc]) , 

eye 

(14) 

where Xfj is constructed with an average pion mass, = ^m n o + ^m n ± , and J2 C yc 
is a sum over the three cyclic exchanges of nucleons i,j,k. For the Urbana models 
Q^r = i^2tt> as i n the original Fujita- Miyazawa model fl27|), while other potentials 
like the Tucson-Melbourne (|2^) and Brazil ( ]29| ) models, have a ratio slightly larger 
than i. The shorter-range phenomenological term is given by 

K?* = E Afl^CTivryjr 2 ^^) . (15) 

eye 

For the Urbana IX (UIX) model, the parameters were determined by fitting the 
binding energy of 3 H and the density of nuclear matter in conjunction with AV18. 

As shown below, while the combined AV18/UIX Hamiltonian reproduces the 
binding energies of s-shell nuclei, it does not do so well for light p-shell nuclei. 
Recently a new class of NNN potentials, called the Illinois models, has been de- 
veloped to address this problem (p2;). These potentials contain the Urbana terms 
and two additional terms, resulting in a total of four coupling constants that can 
be adjusted to fit the data. 

One term, V^ 1 , is due to irN S-wave scattering as illustrated in Fig. |]b. It 
has been included in a number of NNN potentials like the Tucson-Melbourne (p8p 
and Brazil ( p9|) models. The Illinois models use the form recommended in the 
latest Texas model (30), where chiral symmetry is used to constrain the structure 



of the interaction. However, in practice, this term is much smaller than the 
contribution and behaves similarly in light nuclei, so it is difficult to establish its 
strength independently just from calculations of energy levels. 

A more important addition is a simplified form for three-pion rings containing 
one or two deltas (Fig. |l]c,d). As discussed in Ref. (|22|), these diagrams result in 
a large number of terms, the most important of which are used to construct the 
Illinois models: 

] 50 ql ql , 26 Al aI 

3 iJ T O a -f 3 -n-T-ft-n 



t/ 3tt,A/? _ A AR 
V ijk — Ji ^r 



(16) 



Here the S^. and A^. are operators that are symmetric or antisymmetric under any 
exchange of the three nucleons, and the subscript x = o or r indicates that the 
operators act on, respectively, spin and space or just isospin degrees of freedom. 
The S\ is a projector onto isospin-| triples: 

Si = 2 + f (Ti ■ Tj + Tj -T k + T k - Ti) = 4P T=3/2 . (17) 

To the extent isospin is conserved, there are no such triples in the s-shell nuclei, 
and so this term does not affect them. It is also zero for Nd scattering. However, 
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the S^Sp term is attractive in all the p-shell nuclei studied. The A\. has the same 
structure as the isospin part of anticommutator part of V 2lT ' P , but the A^A^ term 
is repulsive in all nuclei studied so far. In p-shell nuclei, the magnitude of the 
A^A^ term is smaller than that of the S^S^ term, so the net effect of the V^,' AR is 
slight repulsion in s-shell nuclei and larger attraction in p-shell nuclei. The reader 
is referred to Ref. (^), and its appendix, for the complete structure of V^' R . 
Relativistic corrections to the Hamiltonian of Eq.(|l]) have been studied in three- 



and four-body nuclei (31,32,B3). We only briefly outline the results here. If a 



square-root kinetic energy, 



Kr* = Jti + m?-mi, (18) 



where pi is the momentum of nucleon i, is introduced into the Hamiltonian, then 
the AW potential must be readjusted to refit the NN data. This can be done for 
AV18 with only small adjustments of the parameters, to produce a v rel . VMC 
calculations show that the resulting (K rel + v rel ) ~ (K + v) for 3 H and 4 He; thus 
the square-root kinetic energy may be neglected in light nuclei. 

A second relativistic correction, the so-called "boost correction," arises from the 
fact that the two-nucleon interaction Vij depends both on the relative momentum 
p = (pi — Pj)/2 and the total momentum P = pj + pj of the interacting nucleons, 

%-(p,P) =%-(p) +Sv ij {p,P) , (19) 

where 5vij(p, P = 0) = 0. Fitting the two-body potential to two-nucleon scat- 
tering data determines only v^j. The 5vij can be determined from Vif, a suitable 
approximation (33,34,35|) is 



where only the first six terms (the static terms) of retained and the gradient 

operators do not act on Vij. The 6vij is meant to be used only in first-order 
perturbation in a non-relativistic wave function. Its expectation values in light 
nuclei are nearly proportional to (V^ k ) (p2[) , and thus one can consider that part 
of Aji, Eq. (|T5|) , comes from 5vij. This proportionality does not, however, extend 
to nuclear matter or pure neutron systems. There one must use a reduced A^ that 
does not contain the part ascribable to 5vij, and explicitly evaluate (5%); doing 
so results in a softer equation of state for dense matter. 



3 VARIATIONAL MONTE CARLO 



The variational method can be used to obtain approximate solutions to the many- 
body Schrodinger equation, = Efy, for a wide range of nuclear systems, 
including few-body nuclei, light closed-shell nuclei, nuclear matter, and neutron 
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stars (36). A suitably parameterized wave function, tyy, is used to calculate an 
upper bound to the exact ground-state energy, 

(V V \H\9 V ) 



Ey 



> E n 



(21) 



(*v|*v) 

The parameters in tyy are varied to minimize Ey, and the lowest value is taken as 
the approximate solution. Upper bounds to excited states are also obtainable, if 
they have different quantum numbers from the ground state, or from small-basis 
diagonalizations if they have the same quantum numbers. The corresponding *Sfy 
can be used to calculate other properties, such as particle density or electromagnetic 
moments, or to start a Green's function Monte Carlo calculation. In this section we 
describe the ansatz for tyy for light nuclei and briefly review how the expectation 
value is evaluated and the parameters of ^y are fixed. 
The best variational wave function has the form (pd|) 



\*v) 



1 + E (U ijh + U™ 1 ) + 



i<j<k 

where the pair wave function, ^p, is given by 



i<3 



LS 



l*P> , 



S]J(1 + U, 



I.I' 



I* 



(22) 



(23) 



The Uij, U(j S , Uijk, and U™ 1 are noncommuting two- and three- nucleon correla- 
tion operators, and S is a symmetrization operator. The form of the antisymmetric 
Jastrow wave function, tyj, depends on the nuclear state under investigation. For 
the s-shell nuclei the simple form 



II fijk( r ik,r jk )Y[f c (r 



\$ A (JMTT 3 )) 



(24) 



i<j<k i<j 

is used. Here f c (i~ij) and ff^ k are central (spin-isospin independent) two- and three- 
body correlation functions and <3?a is an antisymmetrized spin-isospin state, e.g., 



\^(\\\\))=A\}p[p}n) , 
|$ 4 (0000)) =A\\pip\n{n) 
with A the antisymmetrization operator. 



(25) 
(26) 



The two-body correlation operators Uij and U^ s are sums of spin, isospin, tensor, 



and spin-orbit terms: 



U 



U, 



i.i 



LS 
ij 



E 

p=2,6 

E 

p=7,8 



IT fijk( T ik: r jk) 



u p (rij)0: 



n fijki^k^jk) 



k^ij 



(27) 
(28) 
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Figure 2: Correlation functions in light nuclei: the central pair correlation, f c , 
and tensor-isospin correlation, ut T , for s-shell nuclei are shown with solid lines. 
The dashed lines show /^( 3 H) and /^( 4 He) to illustrate the long-range separation 
behavior. The dot-dashed lines show the central correlation between s- and p-shell 

LS\n] r 

nucleons, f sp , the ut T , and the one-body wave function, <p p , for Li. 

where the Ofj were introduced in Eq.(|Tl|). The spin-orbit correlations are only 
summed in Eq.(|22|) because of the extra computational expense of evaluating pow- 
ers of L • S that would occur if it was inserted in the symmetrized product of 
Eq.(H). 

The central / c (r) and noncentral u p (r) pair correlation functions reflect the in- 
fluence of the two-body potential at short distances, while satisfying asymptotic 
boundary conditions of cluster separability. Reasonable functions are generated by 
minimizing the two-body cluster energy of a somewhat modified NN interaction; 
this results in a set of eight coupled, Schrodinger-like, differential equations corre- 
sponding to linear combinations of the first eight operators i of Vij, with a number 
of embedded variational parameters (|To|). The / c (r) is small at short distances, to 
reduce the contribution of the repulsive core of v^, and peaks at an intermediate 
distance corresponding to the maximum attraction of v^, as illustrated in Fig. ^ 
for several nuclei. For the s-shell nuclei, f c (f) falls off at larger distances to keep 
the system confined. For example, in He, [/ c (^)] 3 ~ exp(— nr)/r at large r, as 
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shown by the dashed line in Fig. ^, where k corresponds to an a — > t-\-p separation 
energy. The noncentral u p (r) are all relatively small; the most important is the 
long-range tensor-isospin part ut T (r), also shown in Fig. ^ for several nuclei, which 
is mainly induced by the one-pion-exchange part of Vij. 

The ffj k , and Uijk are three-nucleon correlations also induced by Vij (|TT|). 
The U™ 1 are three-body correlations induced by the three-nucleon interaction, 
which have the form suggested by perturbation theory: 

U ijk* = E e x V ijk{fij,f jk , hi) , (29) 

X 

with f = yr, y a scaling parameter, and e x a (small negative) strength parameter. 
A somewhat simpler wave function to evaluate than ^>y is given by: 

|* T ) = 

where UEk is a slightly truncated TNI correlation. The gives about 1 — 2 
MeV less binding than the full but costs less than half as much to construct, 
making it a more efficient starting point for the GFMC calculation fll8|). 

The Jastrow wave function for the light p-shell nuclei is more complicated, as 
a number of nucleons must be placed in the unfilled p-shell. The LS coupling 
scheme is used to obtain the desired JM value of a given state, as suggested by the 
shell-model studies of light p-shell nuclei (|5|). Different possible LS combinations 
lead to multiple components in the Jastrow wave function. The possibility that 
the central correlations / c (j"y) could depend upon the shells (s or p) occupied by 
the particles and on the LS coupling is also allowed for: 

= a\ n ftjk n fssirij) n /*(*•«) (31) 

[i<j<fc *<i<4 k<4<l<A 

E {Plsw II fp P S[n \ri m )\^A(LS[n]JMTT 3 ) 1234:56 ... A )) 1 . 

LS[n] A<l<m<A J 

The operator A indicates an antisymmetric sum over all possible partitions of the 
A particles into 4 s-shell and (A — 4) p-shell ones. The central correlation f ss (r) is 
the f c (f) from the 4 He wave function. The f sp (r), shown in Fig. [2| for Li nuclei, is 
similar to the / c (r) at short range, but with a long-range tail going to unity; this 
helps the wave function factorize to a cluster structure like a + d in 6 Li or a + 1 in 
7 Li at large cluster separations. The fpp (r) is similar to the deuteron (triton) 
f c (r) in the case of 6 Li ( 7 Li). 

The LS components of the single-particle wave function are given by: 

\$A(LS[n]JMTT 3 ) 12 u:56...A) = |$a(0000) 1234 } II ^ 5W (^) (32) 

4<1<A 



1 + 



e 

i<j<k 



I*/ 



(30) 
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Table 1: Permutation symmetry terms for LS-coupling in A = 6 — 8 nuclei and 
corresponding spin states. 



A [re] L (T, S) highest symmetry states 



6 [2] 


o, 


2 


(1,0)(0,1) 


t, He(0+,2+), D Li(l+ 2+ 3+) 


[11] 


1 




(1,1)(0,0) 


6 He(l+) 


7 [3] 


1, 


3 


V2' 2> 


1J1 \2 '2 '2 '2 ' 


[21] 


1, 


2 


(|> ' 2/ 
(f' 2^2 ' 2~) 




[111] 









8 [4] 


0, 


2,4 


(0,0) 


8 Be(0+,2+,4+) 


[31] 


1. 


2,3 


(1,1)(1,0)(0,1) 


8 Li(0+ -4+), 8 Be(l+,3+) 
8 He(0+,2+) 


[22] 


o. 


2 


(2,0)(1,1)(0,2)(0,0) 


[211] 


1 




(2,1)(1,2)(1,1)(1,0)(0,1) 


8 He(l+) 



\[ II ^((^aOllMiln] X [ JJ Xz(|™*)]SM s > ><[ ^(^ 3 )]tT 3 > • 

[ 4<«<A 4</<A J JM 4<«<A 

The (t>p S ^ n \R a k) are p-wave solutions of a particle in an effective a + N potential 
and are functions of the distance between the center of mass of the a core and 
nucleon k; they may be different for different L5[n] components. A Woods-Saxon 
potential well, 

V aN (r) = V p LS [l + expC-^T 1 , (33) 

v a p 

where V p LS , R p , and a p are variational parameters, and a Coulomb term if appro- 
priate is used. 

For A > 7 nuclei, with three or more p-shell particles there are multiple ways of 
coupling the orbital angular momentum, L, spin, S, and isospin, T, to obtain the 
total (J n ;T) of the system. The spatial permutation symmetry, denoted by the 
Young pattern [n], is used to enumerate the different possible terms, as discussed 
in Appendix 1C of Ref. (|37|) . The different possible contributions to A = 6-8 
nuclei are given in Table [y, with the corresponding spin states of highest spatial 
symmetry for each nucleus. 

After other parameters in the trial function have been optimized, a series of 
calculations are made in which the @LS[n] °1 Eq.(p2|) may be different in the left- 
and right-hand-side wave functions to obtain the diagonal and off-diagonal matrix 
elements of the Hamiltonian and the corresponding normalizations and overlaps. 
The resulting NxN matrices are diagonalized to find the Pis[n] eigenvectors, using 
generalized eigenvalue routines because the correlated ^>v are n °t orthogonal. This 
allows us to project out not only the ground states, but excited states of the same 
(J 77 ; T) quantum numbers. For example, the ground states of 6 Li, 'Li, and 8 Li are 
obtained from 3x3, 5x5, and 7x7 diagonalizations, respectively. 



Exact Calculations of Light Nuclei 



13 



The energy expectation value of Eq.(|2l|) is evaluated using Monte Carlo inte- 
gration. A detailed technical description of the methods used can be found in 



Refs. (fLC| 38). Monte Carlo sampling is done both in configuration space and 
in the order of operators in the symmetrized product of Eq.(|23|) by following a 
Metropolis random walk. The expectation value for an operator O is computed 
with the expression 



(O) 



Z P>q JdR~ 


*t(R) * g (R)/W p? (R) 


' W pq (K) 


J2p,q I ^R 


"4(R)* ? (R)/W P? (R)" 


W pq (K) 



(34) 



where samples are drawn from a probability distribution, W pq (Fl). The subscripts p 
and q specify the order of operators in the left- and right-hand-side wave functions, 
while the integration runs over the particle coordinates R = (ri,r2, . . . , r^). The 
probability distribution is constructed from the f p of Eq. (p3) : 



(35) 



W pq (R) = |Re((*^ p (R)* Pj9 (R; 

This is much less expensive to compute than using the full wave function of Eq.(p2|) 
with its spin-orbit and operator-dependent three-body correlations, but it typically 
has a norm within 1-2% of the full wave function. 

Expectation values have a statistical error which can be estimated by the stan- 
dard deviation a: 



a 



(Q 2 ) - {Of 

N - 1 



1/2 



(36) 



where ./V is the number of statistically independent samples. Block averaging 
schemes can also be used to estimate the autocorrelation times and determine 
the statistical error. 

complex numbers, 
(37) 



The wave function \& can be represented by an array of 2 J 

*(R)=£> a (R)|a) , 



(I) 



where the V'a(R) are the complex coefficients of each state \a) with specific third 
components of spin and isospin. This gives vectors with 96, 1280, 4480, and 14336 
complex numbers for 4 He, 6 Li, 7 Li, and 8 Li, respectively. The spin, isospin, and 
tensor operators Of~ 2 ' 6 contained in the two-body correlation operator Uij, and in 
the Hamiltonian are sparse matrices in this basis. For forces that are largely charge- 
independent, as is the case here, this charge-conserving basis can be replaced with 
an isospin-conserving basis that has N(A, T) = 2 A x I(A, T) components, where 



I(A,T) 



2T + 1 
\A + T + l 



A 

\A + T 



(38) 



This reduces the number of vector elements to 32, 320, 1792, and 7168 for the 
cases given above — a significant savings. In practice, isospin operators are more 
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expensive to evaluate in this basis, but the overall savings in computation is still 
large. Furthermore, if for even- A nuclei the M = state is computed, only half the 
spin components need to be evaluated; the other ones can be found by time-reversal 
invariance. 

Expectation values of the kinetic energy and spin-orbit potential require the com- 
putation of first derivatives and diagonal second derivatives of the wave function. 
These are obtained by evaluating the wave function at 6A slightly shifted positions 
of the coordinates R and taking finite differences, as discussed in Ref. (|l0|). Poten- 
tial terms quadratic in L require mixed second derivatives, which can be obtained 
by additional wave function evaluations and finite differences. 

The rapid growth in the size of VP as A increases is the chief limitation, both 
in computer time and memory requirements, on extending the current method 
to larger nuclei. The cost of an energy calculation for a given configuration also 
increases as the number of pairs, P = \A{A — 1), because of the number of pair 
operations required to construct and roughly as the number of particles, A, due 
to the number of wave function evaluations required to compute the kinetic energy 
by finite differences. These factors are shown in Table || for a number of nuclei, 
with a final overall product showing the cost relative to that for an 8 Be calculation. 
Some initial A = 9, 10 calculations have been made, but are at the limit of current 
computer resources, so while 12 C should be reached in a few years, it may be the 
practical limit for this approach. Calculations of nuclei like 16 or 40 Ca, will require 
other methods such as auxiliary-field diffusion Monte Carlo, which samples the spin 
and isospin components of the wave function by introducing auxiliary fields, as the 
plain Monte Carlo samples the spatial part. Also shown in the table is the cost for 
several pure neutron problems, such as an eight-body drop (with external confining 
potential) or 14 or 38 neutrons in a box (with periodic boundary conditions) as 
neutron matter simulations. The 8 n drop has been calculated (p9|,f4C|) by VMC and 
GFMC, and initial box simulations have been done for 14 neutrons by GFMC (41) 
and up to 54 neutrons by AFDMC @. 

A major problem arises in minimizing the variational energy for p-shell nuclei 
using the above wave functions: there is no variational minimum that gives reason- 
able rms radii. For example, the variational energy for 6 Li is slightly more bound 
than for 4 He, but is not more bound than for separated 4 He and 2 H nuclei, so the 
wave function is not stable against breakup into a + d subclusters. Consequently, 
the energy can be lowered toward the sum of 4 He and 2 H energies by making the 
wave function more and more diffuse. Such a diffuse wave function would not be 
useful for computing other nuclear properties, or as a starting point for the GFMC 
calculation, so the search for variational parameters is constrained by requiring the 
resulting point proton rms radius, r p , to be close to the experimental values for 6 Li 
and Li ground states. For other A = 6-8 ground states, and all the excited states, 
the trial functions contain minimal changes to the 6 Li and 7 Li wave functions, with 
the added requirement that excited states should not have smaller radii than the 
ground states. The final step is always the diagonalization of the Hamiltonian in 
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Table 2: Approximate scaling of VMC and GFMC calculations with system size. 
The number of pairs is P, the state-vector size is the product of the number of 
spin states and number of isospin states, I(A,T); as is discussed in the text, these 
numbers may be halved for even- A nuclei. The last column is the relative difficulty 
scaled to 8 Be, assuming M = states are computed for even- A nuclei. 





A 


P 


2 A x I(A,T) 


n/n( 8 Be) 


3 H 


3 


3 


8x2 


0.0004 


4 He 


4 


6 


16x2 


0.001 


5 He 


5 


10 


32x5 


0.020 


6 Li 


6 


15 


64x5 


0.036 


7 Li 


7 


21 


128x14 


0.66 


8 Be 


8 


28 


256x14 


1. 


9 Be 


9 


36 


512x42 


17. 


io B 


10 


45 


1024x42 


24. 


12 C 


12 


66 


4096x132 


530. 


16Q 


16 


120 


65536x1430 


2 x 10 5 


40 Ca 


40 


780 


7xl0 21 


3 x 10 20 


8 n 


8 


28 


256x1 


0.071 


14 n 


14 


91 


16384x1 


26. 


38 n 


38 


703 


3xlO n 


9 x 10 9 



the PlsWi] mixing parameters. 

4 GREEN'S FUNCTION MONTE CARLO 

The GFMC method (§,|42|) projects out the exact lowest-energy state, f° r a 
given set of quantum numbers, using = nm r-+oo ex p[ — (H — Eq)t]^t, where 
is an initial trial function. If the maximum r actually used is large enough, 
the eigenvalue Eq is calculated exactly while other expectation values are generally 
calculated neglecting terms of order — ^>t\ 2 and higher (|l^). In contrast, the 
error in the variational energy, Ey, is of order — ^t) 2 ) and other expectation 
values calculated with have errors of order — *S>t\- In the following we 
present a brief overview of modern nuclear GFMC methods; much more detail 
may be found in Refs. (18, 19|). 



We start with the of Eq.(pO) and define the propagated wave function v1/(t) 



where we have introduced a small time step, r = nAr; obviously ^(r = 0) = 
and \I/(t — > oo) = The ^(t) is represented by a vector function of R using 
Eq,(|37|), and the Green's function, G Q/ g(R, R') is a matrix function of R and R' 
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in spin-isospin space, denned as 

G a/3 (R,R') = (R,a|e-^- B °) Ar |R',/3) . 



(40) 



It is calculated with leading errors of order (At) 3 . Omitting the spin-isospin indices 
a, f3 for brevity, ^(R n ,r) is given by 



*(Rn,T) = / G(R n ,R n _i)---G(Ri,Ro)*T(Ro)rfRn-i---rfRidR , 



(41) 



with the integral being evaluated stochastically. 

The short-time propagator should allow as large a time step Ar as possible, 
because the total computational time for propagation is proportional to 1/Ar. 
Earlier calculations (17, p|,|42]) used the propagator obtained from the Feynman 
formulae in which the kinetic and potential energy terms of H are separately ex- 
ponentiated. The main error in this approximation to G Qj/ g comes from terms in 
e -HAr ikying mu itiple Vij, like VijKvij(Ar) 3 , where K is the kinetic energy, which 
can become large when particles i and j are very close due to the large repulsive 
core in v^. This requires a rather small At ~ 0.1 GeV -1 . 

It has been found in studies of bulk helium atoms (^) that including the exact 
two-body propagator allows much larger time steps. This short-time propagator is 



G a/3 (R,R') = G (R,R')(a| 



9o,ij( r ij, r 'ij) 



where 



G (R,R') 



(Rle 



—KAt 



R') 



2vr^Ar 



3.4 



exp 



l\2 



-(R-RQ 

2% 2 L\T/m 



gij is the exact two-body propagator, 

9ij ( r ij j r ij ) = ( r 



; e 



-H ij Ar\ J \ 



(42) 



(43) 



(44) 



and <7o,r/ is the free two-body propagator. All terms containing any number of the 
same and K are treated exactly in this propagator, as we have included the 
imaginary-time equivalent of the full two-body scattering amplitude. It still has er- 
rors of order (Ar) 3 , however they are from commutators of terms like ^K%(At) 3 
which become large only when both pairs ij and ik are close. Because this is a 
rare occurrence, a five times larger time step, Ar ~ 0.5 GeV -1 , can be used (|T8|). 

Including the three-body forces and the Eq in Eq.(f40|), the complete propagator 
is given by 

G aP (R,R>) = e B - AT Go(R I R')exp[-X;(^?fc(R) + ^?fc(R , ))^ : ] 



<a|/ 3 (R)|7><7l 



9ij( r 



\5){5\h(^!)\(i) , (45) 
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where 



(46) 



and Vjj k = Vjjk + Vj^ represents, in general, all non-central Vijk terms. With the 
exponential of V^ k expanded to first order in At, there are additional error terms 
of the form V^ fe V^j-/ fe /(Ar) . However, they have negligible effect because V?j k has 
a magnitude of only a few MeV. 

Quantities of interest are evaluated in terms of a "mixed" expectation value 
between and ^(r): 



Mixed 



(grjQjg(r)) 
(*t|*(t)) 

/ rfP w ^(R n )OG(R n , Rn_i) • • • G(Ri, Ro)^ r (Ro 
/dP„^(R n )G(R n ,R n _i) • ..G(Ri,Ro)*r(Ro) 



(47) 



where P n = Rq, Ri, • • • , R n denotes the 'path', and dP n = dRo<iRi • • • d¥L n with 
the integral over the paths being carried out stochastically. The desired expectation 
values would, of course, have VP(t) on both sides; by writing \&(t) = + 5^(r) 
and neglecting terms of order [5^(r)] 2 , we obtain the approximate expression 

(0(t)) = \^)Mt)) « (0(T))Mixed + [(0(T))Mixed - (0) T ] , (48) 

where {0)t is the variational expectation value. More accurate evaluations of 
(0(t)) are possible (|44]), essentially by measuring the observable at the mid-point 
of the path. However, such estimates require a propagation twice as long as the 
mixed estimate and require separate propagations for every (O) to be evaluated. 
The nuclear calculations published to date use the approximation of Eq. (|4"q ) . 

A special case is the expectation value of the Hamiltonian. The (H(r))wuxed can 
be re-expressed as 



- (*T\e-l H -*WHe-V-*>Wm 
(H(r)) Mixcd - ^ r | e _ (H _ Bb)r/2e _ (ff _ jEb)T/2 |^ T) > , (49) 

since the propagator exp[—(H — Eq)t] commutes with the Hamiltonian. Thus 
(-fT(r))Mixed approaches Eq in the limit r — > oo, and the expectation value obeys 
the variational principle for all r. 

The AVI 8 interaction contains terms that are quadratic in orbital angular mo- 
mentum, L. These terms are, in essence, state- and position-dependent modifica- 
tions of the mass of the nucleons. If they are included in the calculation of gij, 
then the ratio gij/go^ij in Eq.(43) becomes unbounded for large |ry| or |r£ -| 



Monte Carlo statistical error will also be unbounded. Hence the GFMC propaga- 
tor is constructed with a simpler isoscalar interaction, H', with a vL that has only 
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eight operator terms, [1, <Tj • crj, Sy, L • S] <g) [1, Tj • Tj], chosen such that it equals 
the CI part of the full interaction in all S- and P-waves and in the deuteron. The 
v'ij is a little more attractive than Vij, so a V^ k adjusted so that (H 1 ) « (H) is also 
used. This ensures the GFMC algorithm will not propagate to excessively large 
densities due to overbinding. Consequently, the upper bound property applies to 
(H'(t)), and (H — H') must be evaluated perturbatively. 

Another complication that arises in the GFMC algorithm is the "fermion sign 
problem". This arises from the stochastic evaluation of the matrix elements in 



Eq.(47). The G(Rj, Rj__i) is a local operator and can mix in the boson solution. 
This has a (much) lower energy than the fermion solution and thus is exponentially 
amplified in subsequent propagations. In the final integration with the antisym- 
metric *&Ti the desired fermionic part is projected out, but in the presence of large 
statistical errors that grow exponentially with r. Because the number of pairs that 
can be exchanged grows with A, the sign problem also grows exponentially with 
increasing A. For ^4>8, the errors grow so fast that convergence in r cannot be 
achieved. 

For simple scalar wave functions, the fermion sign problem can be controlled by 
not allowing the propagation to move across a node of the wave function. Such 
"fixed-node" GFMC provides an approximate solution which is the best possible 
variational wave function with the same nodal structure as VPy. However, a more 
complicated solution is necessary for the spin- and isospin-dependent wave func- 
tions of nuclei. In the last few years a suitable "constrained path" approximation 
has been developed and extensively tested, first for condensed matter systems ( fl6| ) 
and more recently for nuclei fllS|). The basic idea of the constrained-path method 
is to discard those configurations that, in future generations, will contribute only 
noise to expectation values. If the exact ground state |^o) were known, any con- 
figuration for which 

tt(P„)ttf (Rn) = , (50) 

where a sum over spin-isospin states is implied, could be discarded. Here P n 
designates the complete path, (Rq, ...,R n ) that has led to R n . The sum of these 
discarded configurations can be written as a state Iv&d), which obviously has zero 
overlap with the ground state. The contains only excited states and should 
decay away as r — > oo, thus discarding it is justified. Of course the exact is 
not known, and so configurations are discarded with a probability such that the 
average overlap with the trial wave function, 

(* d |* T )=0. (51) 



Many tests of this procedure have been made fll9| ) and it usually gives results that 
are consistent with unconstrained propagation, within statistical errors. However 
a few cases in which the constrained propagation converges to the wrong energy 
(either above or below the correct energy) have been found. Therefore a small 
number, n u = 10 to 20, unconstrained steps are made before evaluating expectation 
values. These few unconstrained steps, out of typically 400 total steps, appear to 
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Figure 3: GFMC E{r) using AV18/UIX for various states of 8 Be. Averages with 
their standard errors are shown by the solid and dashed lines. 



be enough to damp out errors introduced by the constraint, but do not greatly 
increase the statistical error. 

Figure || shows the progress with increasing r of typical constrained GFMC 
calculations, in this case for various states of 8 Be. The values shown at r = 
are the VMC values using ^>t] VMC values with the best known ^fy are about 
1.5 MeV below these. The GFMC very rapidly makes a large improvement on 
these energies; by r = 0.01 MeV -1 , the energies have been reduced by 6.5 to 
8 MeV. This rapid improvement corresponds to the removal of small admixtures of 
states with excitation energies ~ 1 GeV from Averages over typically the last 
nine r values are used as the GFMC energy. The standard deviation, computed 
using block averaging, of all of the individual energies for these r values is used to 
compute the corresponding statistical error. The solid lines show these averages; 
the corresponding dashed lines show the statistical errors. The g.s., 1 + , and 3 + 
states of 8 Be have widths less than 250 keV and their E(r) appear to be converged 
and constant over the averaging regime. The 2 + state has a width of 1.5 MeV and 
the 4 + state's width is 3.5 MeV; the E(t) for the 2 + state might be converged, but 
that of the 4 + state is clearly steadily decreasing. Reliable estimates of the energies 
of broad resonances requires the use of scattering-state boundary conditions (see 
Sec. 0) which have not yet been implemented for more than five nucleons. 

Tests with different as starting points for the GFMC calculation are described 
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in Ref. (18). The most crucial aspect in choosing is that it have the correct 
mix of spatial symmetries, i.e., the PlsWi] admixtures. This is because the limited 
propagation time is not sufficient for the GFMC algorithm to filter out low-lying 
excitations with the same quantum numbers. This issue is unimportant in 4 He, 
where the first + excited state is near 20 MeV, but in 6 Li, the first 1+ excited 
state is at only 5.65 MeV; other p-shell nuclei have similar low-lying excitations. 
Otherwise, the GFMC algorithm is able to rapidly correct for some very poor 
whose variational energies are actually positive. However, a good helps to 
keep the error bars small at larger r. The most efficient balance between speed 
of construction of and smallest number of samples needed to achieve a given 
error bar seems to be given by Eq.(|30|). A number of other tests of the GFMC 



algorithm and its ability to determine nuclear radii are described in Refs. ( [1S| , 19). 

The calculations described here are computationally intense, and would not have 
been possible without the advent of parallel supercomputers. The initial stud- 



ies (17) of Li required ~ 2,000 node hours on an IBM SP1 to propagate 10,000 
configurations to r = 0.06 MeV" 1 . Today, with improvements in the algorithm 
and after much effort to optimize the computer codes, the same calculation requires 
about 15 node hours on a third generation IBM SP. However, a present calculation 
for 8 Li, requiring 10,000 configurations to get a reasonable error bar, takes 1,000 
node hours, so forefront computer resources remain essential for this program. 

5 ENERGY RESULTS 

There are a number of accurate many-body methods for evaluating the bind- 
ing energy of the s-shell nuclei, 3 H, 3 He, and 4 He, using realistic nuclear forces. 
These include Faddeev in configuration space (^7|) (FadR), Faddeev and Faddeev- 
Yakubovsky in momentum space (^) (FadQ), and hyperspherical harmonics ( |49| ) 
(HH), correlated hyperspherical harmonics (|50| ) (CHH) and pair-correlated hyper- 
spherical harmonics (^Tj) (PHH). Some results of these methods for the AV18 and 
AV18/UIX Hamiltonians are shown in Table ||. We observe that the VMC upper 
bounds are generally 1.5-2.0% above the GFMC energies, which in turn are 0.25- 
1.0% above the Faddeev results. The disagreement between GFMC and the other 
methods may be attributable to the fact that GFMC is propagated with an H' and 
a small piece of the Hamiltonian is computed in perturbation. At this fine level of 
comparison, one needs to worry about how features such as T = | admixtures and 
n—p mass differences in the trinucleon ground state are treated in each calculation. 

Table || also shows the necessity of including (at least) a three-nucleon potential 
in the Hamiltonian in order to reproduce the A=3,4 experimental binding energies. 



This is true for all the other modern AW potentials, as shown in Refs. (47, 48). 
The local potentials AV18, Reid93, and Nijm II, give very similar results, while the 
slightly nonlocal Nijm I gives 1-3% more binding, and the more nonlocal CD-Bonn 
gives 5-8% more binding, or about halfway between the local potential values and 
experiment. At present, no two- plus three-nucleon potential combination gives 
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Table 3: Energies of 3 H and 4 He in MeV, computed by VMC and GFMC compared 
to momentum-space Faddeev (fS|), configuration-space Faddeev ( [47|) and various 
hyper spherical harmonics methods (^). The HH column values are PHH for 3 H; 
HH for 4 He, AV18; and CHH for 4 He, AV18/UIX 



Hamiltonian 


nucleus 


VMC 


GFMC 


FadQ 


FadR 


HH 


Expt. 


AV18 


3 H 


-7.50(1) 


-7.61(1) 


-7.623 


-7.62 


-7.618 


8.482 




4 He 


-23.72(3) 


-24.07(4) 


-24.28 




-24.18 


28.296 


AV18/UIX 


3 H 


-8.32(1) 


-8.46(1) 


-8.478 




-8.475 


8.482 




4 He 


-27.78(3) 


-28.33(2) 


-28.50 




-28.1 


28.296 



an exact fit to both trinucleon and 4 He energies, but several combinations, like 
AV18/UIX, come quite close. It also appears that it may be easier to fit 3 He and 
4 He simultaneously, rather than 3 H and 4 He, perhaps because of the better data 
constraints on the pp interaction compared to the nn interaction (f48|). In principle, 
there should also be four-body forces, but their contribution must be quite small, 
and it is impossible to unambiguously identify a need for such terms until a more 
thorough survey of possible three-nucleon potentials is made and the discrepancies 
between the various many-body methods are resolved. 

Figure H compares VMC (using the full \Ev of Eq.(|2|)) and GFMC calculations 
for the AV18/UIX Hamiltonian, and also shows experimental energies. All GFMC 
energies in this review are from Ref. (|22j); VMC energies are from Refs. (19, 18). 
All particle-stable or narrow-width (r < 150 keV) states for 4 < A < 8 are shown, 
except isobaric analogs. We see that the VMC calculations get progressively worse 
in the p-shell compared to the s-shell, being on the order of 10-15% above the final 
GFMC results, for nuclei with N ~ Z. In absolute terms, misses roughly an 
extra 1.5 MeV of binding for each p-shell nucleon that is added. The VMC results 
also fail to reproduce important qualitative features of the GFMC calculations; 
for example 6 Li is stable against breakup into a+d for this Hamiltonian but the 
VMC calculation shows it unbound by ~2 MeV (the black dashed lines show the 
indicated thresholds computed using the sub-cluster energies appropriate to each 
calculation or experiment). Clearly, there is some significant feature of light p-shell 
nuclei that is not yet included in the trial wave functions. 

Table shows the GFMC energy values for several Hamiltonians. The table 
and Fig. ^ show that the AV18/UIX Hamiltonian, which gives excellent energies, 
compared to experiment, for the s-shell nuclei, significantly underbinds the p-shell 
nuclei. The Li isotopes are stable against breakup into subclusters with AV18/UIX, 
but progressively more underbound as A increases. There is also an isospin prob- 
lem, in that the He isotopes are even further off from the experimental values, and 
in fact do not show stability against breakup into subclusters. 

The new Illinois three-nucleon potentials were constructed to solve the p-shell 
binding problems. Figure [| and Table || compare GFMC calculations for AV18 
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Figure 4: VMC and GFMC energies using AV18/UIX compared to experiment. 
Black dashed lines show the indicated breakup thresholds for each method. The 
Monte Carlo statistical errors are shown by the light blue and yellow bands. 

with no Vijk and AV18/IL2 to the same experimental energies shown in Fig. |j. 
The AV18/IL2 Hamiltonian does a good job of reproducing these energies; the rms 
deviation from experiment for these levels is only 360 keV, while it is 2.3 MeV for 
AV18/UIX. The AV18 values with no Vijk show the large contribution that three- 
nucleon potentials make to these binding energies; for A = 8 the IL2 increases the 
binding energy by more than 10 MeV. 

Table || shows some splittings for states that in a simple shell-model picture have 
the same L and S but different total J. In all cases the lowest state for each J 
is used. Some of the states, particularly the 5 He states, have large experimental 
widths; calculations such as those described in Sec. [7| would be more reliable. The 
splittings are computed for the AV18, AV18/UIX, and AV18/IL2 Hamiltonians 
and are compared to experimental values. The AV18 with no Vijk significantly 
underpredicts all the splittings except the 5 — I doublet in 7 Li; however in 
this case the relatively large statistical error for the small splitting makes any 
conclusion difficult (the computed splittings are the difference of two independent 
GFMC calculations whose statistical errors must be added in quadrature). The 
AV18/UIX Hamiltonian substantially improves the 5 He splitting, and a VMC study 
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Table 4: GFMC and experimental energies of A=3-8 ground states in MeV. Ener- 
gies of artificially confined 7- and 8-neutron drops are also shown. The experimental 
values are from Refs. ( |S^|54| ,|55| ) , 



AV18 AV18/UIX AV18/IL2 Expt. 



4 He(0+) 


-7 61(1) 


-8.46(1) 


-8.43(1) 


-8.48 


-24.07(4) 


-28.33(2) 


-28.37(3) 


-28.30 


6 He(0+) 


-23.9(1) 


-28.1(1) 


-29.4(1) 


-29.27 


6 Li(l+) 


-26.9(1) 


-31.1(1) 


-32.3(1) 


-31.99 


7 He(|") 


-21.2(2) 


-25.8(2) 


-29.2(3) 


-28.82 


7 Li(l") 
8 He(0+) 


-31.6(1) 


-37.8(1) 


-39.6(2) 


-39.24 


-21.6(2) 


-27.2(2) 


-31.3(3) 


-31.41 


8 Li(2+) 


-31.8(3) 


-38.0(2) 


-42.2(2) 


-41.28 


8 Be(0+) 


-45.6(3) 


-54.4(2) 


-56.6(4) 


-56.50 


8 n(0+) 


-33.47(5) 


-33.2(1) 


-35.8(2) 




-39.21(8) 


-37.8(1) 


-41.1(3) 






Figure 5: GFMC energies using AV18 and AV18/IL2 compared to experiment. 
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Table 5: GFMC and experimental (|55|,|5 
S but different J, in MeV. 



splittings of states with the same L and 









L 


S 


AV18 


AV18/UIX 


AV18/IL2 


Expt. 


5 He 


l- 

2 


3- 
2 


1 


i 

2 


0.6(1) 


1.1(2) 


1.3(2) 


1.20 


6 Li 


2+ 


-3+ 


2 


1 


0.8(1) 


0.9(1) 


2.2(2) 


2.12 


7 Li 


l- 

2 


3~ 
2 


1 


1 

2 


0.5(2) 


0.3(2) 


0.6(3) 


0.47 


7 Li 


5~ 
2 


7~ 
2 


3 


1 
2 


0.7(2) 


0.8(2) 


2.2(3) 


2.05 


8 Li 


1+ 


-2+ 


1 


1 


0.2(4) 


0.6(2) 


1.7(4) 


0.98 


7 n 


3~ 
2 


l- 

2 


1 


1 
2 


1.65(7) 


1.5(1) 


2.8(3) 





showed that an earlier Urbana Vyk significantly increased the computed splitting in 
15 N, resulting in much better agreement with experiment (^6|). However the other 
splittings in the table are still underpredicted with AV18/UIX. The AV18/IL2 
results in good predictions of all the splittings except for the 1 + — 2 + doublet in 
8 Li which is over predicted. 

A detailed breakdown of the GFMC ground-state energies for AV18/IL2 into 
kinetic, two- and three- nucleon interactions is given in Table |(| Because of the 
extrapolation of the mixed expectation values, Eq. (|48l) , these components are not 
as accurate as the total energy, and do not add up to the full amount; indeed 
the sum of the individual kinetic and potential energies differs from (H) by the 
same amount that the GFMC has improved the *$>t energy. Nevertheless, they 
give a good idea about the relative size of the different terms involved. There is a 
big cancellation between the kinetic and two-body potential terms. Consequently, 
while Vijk is less than 8% of Vij in magnitude, its expectation value is up to 50% 
of the net binding (however the net effect of the IL2 Vyk, defined as the difference 
of the binding energies computed without and with V^&, is at most 30%, as can be 
seen in Table |j). This difference between expectation value and net effect is due 
to the large change that induces in (K + vy). 

Among the subcomponents of vy, the one-pion exchange dominates, being 70- 
80% of the total for A > 3 nuclei. Similarly, the two-pion exchange is the 
dominant component of Vijk- The three-pion exchange term is small and repulsive 
in s-shell nuclei, but attractive in the p-shell nuclei, which have T = | triples. 
It is this term which results in the large improvement of the Illinois models over 
the Urbana models for p-shell nuclei; however its contributions are always less 
than 15% of the two-pion Vijk- Finally, we note that the electromagnetic vjj is 
dominated by the Coulomb interaction between protons, but about 17% (8%) of 
its total contribution comes from the magnetic moment and other terms in He (Li) 
isotopes. 

Energy differences among isobaric analog states are probes of the charge- indepen- 
dence-breaking parts of the Hamiltonian. The understanding of these "Nolen- 
Schiffer energies" has been a theoretical problem for 30 years (|58|) . The energies 
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Table 6: Kinetic, K, and potential energy contributions in MeV. Shown are the 
total NN and NNN potential energies, and Vijk] the electromagnetic, vjj, and 
one-pion, t&, parts of Uy-; and the two- and three-pion parts of Vijk- 





K 


Vij 


Vijk 






V'ijk 


^ijk 




51. 


-59. 


-1.5 


0.04 


-45. 


-3.0 


0.18(1) 


4 He 


115.(1) 


-136.(1) 


-8.4(1) 


0.86 


-105. 


-16.3(1) 


0.63(1) 


6 He 


147.(2) 


-171.(2) 


-11.5(3) 


0.87 


-127.(1) 


-20.3(4) 


-0.91(6) 


6 Li 


160.(2) 


-187.(2) 


-11.1(3) 


1.73 


-150.(1) 


-19.8(4) 


-0.44(5) 


7 He 


175.(3) 


-199.(3) 


-16.3(4) 


0.89 


-145.(2) 


-25.(1) 


-3.1(1) 


7 Li 


199.(3) 


-232.(3) 


-14.5(4) 


1.80 


-178.(2) 


-25.6(6) 


-1.1(1) 


8 He 


190.(3) 


-218.(3) 


-16.3(5) 


0.89 


-153.(1) 


-25.6(6) 


-4.0(1) 


8 Li 


242.(2) 


-278.(2) 


-20.6(4) 


1.93 


-211.(1) 


-34.2(5) 


-3.8(1) 


8 Be 


256.(4) 


-303.(3) 


-21.(1) 


3.32 


-234.(2) 


-38.5(9) 


-0.9(2) 


''n 


105.(1) 


-59.(1) 


-3.6(3) 


0.07 


-10. 


-0.1(1) 


-5.4(3) 


8 n 


122.(1) 


-73.(1) 


-3.0(3) 


0.09 


-12. 


0.3(1) 


-5.9(4) 



for a given isomultiplet of states can be expanded as 

E A ,t(T z ) = Y, a n (A,T)Q n (T,T z ) , (52) 

n<2T 

where Q$ = 1, Q\ = T z , and Q2 = |(3T^ — T 2 ) are isoscalar, isovector, and 
isotensor terms (|59[). The isovector and isotensor coefficients a n (A,T), and various 



Table 7: GFMC isovector and isotensor energy coefficients, a n (A,T), [Eq. (|52|)1 
computed using AV18/IL2, in keV. The contributions shown are kinetic energy, 
K CSB ; Coulomb potential, vciipp); other electromagnetic, v 1,R ; and nuclear po- 
tential, v CSB + v CD . The experimental value for 02(8, 1) is extracted using 16.80 
MeV as the excitation of the pure T 



1 isobaric analog state in 8 Be (60). 



a n (A,T) 


K CSB 


vex (pp) 




yCSB _|_ yCD 


Total 


Expt. 


01(3,5) 


14 


649(1) 


29 


64(0) 


757(1) 


764 


oi(6,l) 


16 


1091(5) 


18 


47(1) 


1172(6) 


1173 


02(6,1) 




166(1) 


19 


107(13) 


293(13) 


223 


oi(7 ,|) 


22 


1447(6) 


40 


79(2) 


1588(7) 


1644 


oi(7,|) 


18 


1337(6) 


12 


52(1) 


1420(8) 


1373 


02(7, |) 




137(1) 


7 


36(6) 


180(7) 


175 


oi(8,l) 


23 


1686(5) 


24 


76(1) 


1810(6) 


1770 


02(8,1) 




141(1) 


4 


-3(8) 


143(8) 


128 


ai(8,2) 


18 


1528(7) 


17 


59(1) 


1622(8) 


1659 


02(8,2) 




136(1) 


6 


38(5) 


180(5) 


153 



2(3 
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Table 8: VMC isospin-mixing matrix elements for 8 Be in keV computed using 



J* 


K OSB 


vciipp) 






Eqi 


Expt. 


2+ 


2 


62 


19 


26 


109(4) 


149 


1+ 


1 


39 





15 


55(2) 


120 


3+ 


1 


33 


15 


13 


62(2) 


63 



contributions to them, are given in Table |. These were calculated as expecta- 
tion values in the T z = —T GFMC wave functions for AV18/IL2. The individual 
terms are: 1) the effect of the neutron-proton mass difference on the kinetic en- 
ergy, K CSB ] 2) the proton-proton Coulomb potential, including the AV18 form 
factor, vciipp); 3) all other electromagnetic terms such as vacuum polarization 
and magnetic-moment terms, -u 7 '^; and 4) the strong-interaction contributions, 
V CSB w ] 1 j c i 1 contributes to the isotensor coefficients and v which contributes to 
the isovector coefficients. 

The 3 H- 3 He mass difference is 757 keV with the AV18/IL2 Hamiltonian, in excel- 
lent agreement with the experimental value of 764 keV. The bulk of the difference 
is the pp Coulomb energy, but 108 keV comes from the other terms. Because 
of its very long range, the vci(pp) dominates the a\{A,T) of heavier nuclei even 
more and small errors in the rms radius of the nucleus can result in changes in the 
vciipp) contribution that can totally mask the effects of the other, much smaller, 
terms. 

Two (2 + ;0+l) states occur very close together in the spectrum of 8 Be at 16.6 
and 16.9 MeV excitation; these isospin-mixed states come from blending the (2 + ;l) 
isobaric analog of the 8 Li ground state with the second (2 + ;0) excited state. There 
are also fairly close (1 + ;0,1) and (3 + ;0,l) pairs at slightly higher energies in the 8 Be 
spectrum. The isospin-mixing matrix elements that connect these pairs of states, 

E m (J) = (tf(J+; 0)|fl-|*(J+; 1)) , (53) 

have been computed using VMC wave functions for the AV18/UIX Hamiltonian. 
Results for the E^iiJ) are given in Table |8[ The experimental values are determined 
from the observed decay widths and energies (p0|). The dominant contribution, 
from the Coulomb potential, typically accounts for less than half of the matrix 
element. We see that the remaining part of the electromagnetic interaction and the 
strong CSB interaction can provide a significant boost, although the experimental 
mixing is still underpredicted by ~20%. It appears that these mixing elements are 
a more sensitive test of the small CSB components of the Hamiltonian than are 
the isobaric analog energy differences. 
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Table 9: GFMC values for point proton rms radii (in fm) and quadrupole moments 
(in fm 2 ) in impulse approximation, computed with AV18/IL2. Experimental values 





(r 2 p ) 


1/2 


Q 






GFMC 


Expt. 


GFMC 


Expt. 


'-Itt 

d R 


1.59(0) 


1.60 






3 He 


1.76(0) 


1.77 






4 He 


1.45(0) 


1.47 






6 He 


1.91(1) 








6 Li 


2.39(1) 


2.43 


-0.32(6) 


-0.083 


7 Li 


2.25(1) 


2.27 


-3.6(1) 


-4.06 


7 Be 


2.44(1) 




-6.1(1) 




8 He 


1.88(1) 








8 Li 


2.09(1) 




3.2(1) 


3.11(5) 


8 B 


2.45(1) 




6.4(1) 


6.8(2) 



6 DENSITY AND MOMENTUM DISTRIBUTIONS 

The one- and two-nucleon distributions of light p-shell nuclei are interesting in a 
variety of experimental settings. For example, the Borromean 6 He and 8 He nuclei 
are popular candidates for study as 'halo' nuclei whose last neutrons are weakly 
bound. In addition, the polarization densities of 6 Li and 7 Li are important because 
of possible applications in polarized targets. 

Point proton rms radii and quadrupole moments are shown in Table |[J, as com- 
puted by GFMC for the AV18/IL2 model, along with experimental values. The 
experimental charge radii have been converted to point proton radii by removing 
the proton and neutron (r 2 ) of 0.743 and -0.116 fm 2 , respectively. The radii are 
in good agreement with experiment, thanks at least in part to the fact that the 
AV18/IL2 model reproduces the binding energies of these nuclei very well. The 
quadrupole moments have been calculated in impulse approximation; two-body 
charge contributions are expected to provide only a few percent correction. The 
agreement with experiment is again fairly good, with the exception of the 6 Li 
quadrupole moment, which involves a delicate cancellation between the contri- 
butions from the deuteron quadrupole moment and the D-wave part of the a+d 
relative wave function. In general, quadrupole moments are difficult to calculate 
accurately with quantum Monte Carlo methods because they are dominated by 
the long-range parts of the wave functions, which contribute very little to the total 
energy that VMC and GFMC are both designed to optimize. 

Calculated and experimental isoscalar and isovector magnetic moments are shown 



in Table 10, as defined with the convention used in Eq.(52); thus the isovector val- 



ues for the T = \ cases are —2 times those often quoted in the literature. (For the 
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Table 10: GFMC values for isoscalar and isovector magnetic moments (in u-n) in 
impulse approximation using AV18/IL2. Experimental values are from Refs. (|53| , 





T 


Isoscalar 


Isovector 






GFMC Expt. 


GFMC Expt. 


3 He- a H 


i 

2 


0.403(0) 0.426 


-4.330(1) -5.107 


6 Li 





0.817(1) 0.822 




7 Be- 7 Li 


1 

2 


0.894(1) 0.929 


-3.93(1) -4.654 


8 B- 8 Li 


1 


1.276(1) 1.345 


0.369(9) -0.309 



A = 8, T = 1 nuclei the experimental isoscalar and isovector moments are obtained 
from the sum and difference of the values for B and Li, since the magnetic moment 
of the T = 1 J 71 " = 2 + state in Be is not measured.) The moments were calculated 
as expectation values in the GFMC wave functions with T z = —T in impulse ap- 
proximation. These impulse isoscalar moments are quite close to experiment and 



it is expected that two-body current contributions are small (64). For isovector 
moments, however, there can be significant corrections from meson-exchange con- 
tributions. We note that the corrections needed for the ^4=7,8 isovector moments 
are the same sign and magnitude as for ^4=3, so an exchange-current model that 
fixes the s-shell ( |65| ) is likely to work for the light p-shell nuclei also. 

In Fig. |6|, we present proton and matter (proton+neutron) densities for the 
stable helium isotopes, as calculated with GFMC for AV18/IL2. Nucleon densities 
are calculated as simple 5-function expectation values, with possible spin and/or 
isospin projectors: for example, the proton density is given by 

P P (r) = 4^2 (^IE 1 -^6(r-\r l -R cm \)\*} . (54) 

i 

The blue symbols and curves show results for 6 He, red ones for 8 He; open symbols 
give the proton distributions which are also interpreted as the alpha "core" density, 
full symbols are the total matter densities. It can be seen that, as more neutrons 
are added, the tails of the matter distributions broaden considerably because of the 
relatively weak binding of the p-shell neutrons. In addition, the central neutron 
and proton densities decrease rather dramatically. This effect does not necessarily 
require any changes to the a core, but can be understood at least partially from 
the fact that the a no longer sits at the center of mass of the entire system. The 
motion relative to the center of mass spreads out the mass distribution relative to 
that of 4 He. 

The figure also shows two attempts to extract the densities of 6,8 He from scatter- 
ing of beams of these short-lived nuclei from proton targets (^6[). The analysis of 
Alkhazov, et al. (dashed curves) used Glauber theory and assumed proton (core) 
and matter density distributions which were varied to fit the measured cross sec- 
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Figure 6: Proton and matter densities in 4 He, 6 He, and 8 He. Analyses of ^He- 
proton scattering are shown by the dashed (Ref. (|66|)) and dot-dashed (Ref. (p?|)) 
curves. 



tions. The later work of Al-Khalili and Tostevin ( |67|) (dot-dashed curves) improved 
on this by doing the Glauber theory using model wave functions which contained 
correlations between the valence neutrons. The GFMC calculations definitely pre- 
fer the latter analysis. 

The proton-proton distribution function is defined by 

M0 = ^(*lE 1 Ap^ l A^ i 5(r-\v i -r j \m . (55) 

These distributions are directly related to the Coulomb sum measured in inclusive 
longitudinal electron scattering; such measurements in 3 He have been used to put 
constraints on the Pppifij), and realistic calculations agree with the experimental 
results(|68[). The behavior of p pp (r) at short distances is largely determined by the 
repulsive core of the NN potential and is nearly independent of the nucleus, but 
at larger distances it is determined by the size of the nucleus. 

It is interesting to compare p pp for 4 He to that of 6 He and 8 He. These distribution 
functions are shown in Fig. [?], again calculated with GFMC for the AV18/IL2 
model. These nuclei each have just one pp pair which presumably is in the "alpha 
core" of 6 ' 8 He. Unlike the one-body densities, these distributions are not sensitive 
to center of mass effects, and thus if the alpha core of 6 ' 8 He is not distorted by 
the surrounding neutrons, all three p pp distributions in the figure should be the 
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Figure 7: Proton-proton distribution functions in helium isotopes 



same. We see that the pp distribution spreads out slightly with neutron number 
in the helium isotopes, with an increase of the pair rms radius of approximately 
4% in going from 4 He to 6 He, and 8% to 8 He. While this could be interpreted as 
a swelling of the alpha core, it might also be due to the charge-exchange (r^ • Tj) 
correlations which can transfer charge from the core to the valence nucleons. Since 
these correlations are rather long-ranged, they can have a significant effect on the 
pp distribution. VMC calculations of 4 He with wave functions modified to give 
p pp distributions close to those of 6 ' 8 He suggest that the alpha cores of 6 ' 8 He are 
excited by ~ 80 and ~ 350 keV, respectively. 

In general, VMC calculations give one-body densities very similar to the GFMC 
results, although the two-body densities may differ by up to 10% in the peak, 
with the GFMC having a somewhat sharper structure. Because the VMC wave 
functions are simpler and less expensive to construct they have been used in a 
number of applications where the single-particle structure is dominant. All these 
calculations were carried out with wave functions for the AV18/UIX Hamiltonian, 
but we expect that GFMC calculations with the improved AV18/IL2 model will 
not qualitatively alter the results. (Remember that the VMC wave functions are 
constrained to reproduce available experimental rms radii.) 

VMC calculations of elastic and inelastic electromagnetic form factors for 6 Li are 
shown in Fig. ||. These have been made in impulse approximation (IA) and with 
meson-exchange contributions (MEC) to the charge and density operators (|l2|). 
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Figure 8: Calculated and experimental elastic and transition form factors in 6 Li. 



Experimental data are from Refs. (W,mL 71, 72, 73, 74, 75). 



The comparison with data for the elastic longitudinal form factor, Fi(q 2 ), is ex- 
cellent, as is the E2 transition to the 3 + first excited state. The MEC corrections 
are small, but stand out noticeably in the first minimum where they significantly 
improve the fit to data. The elastic transverse form factor, Fx{q 2 ), is good up to 
the first zero, but is a little too large in the region of the second maximum. The Ml 
transition to the + isobaric analog state is also reproduced reasonably well. This 
kind of quantitative agreement with data has not been achieved in the past with 
either shell model or ct+d cluster wave functions. In fact, shell model calculations 
normally require the introduction of effective charges, typically adding a charge of 
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Figure 9: Differential cross sections for 7 Li(-7r, n') scattering at = 164 MeV. 
Data are from Ref. (|77|). 



~ 0.5e to both the proton and neutron, to obtain reasonable transition rates (76). 
No effective charges are used in the VMC calculations. 

While electron scattering from nuclei is primarily sensitive to the proton dis- 
tributions in nuclei, ir~ scattering is most sensitive to the neutron distributions. 
Pion inelastic scattering at medium energy (80 < Ei a j, < 300 MeV) is dominated 
by strong absorption due to excitation of the A resonance, and is well described in 
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Figure 10: Experimental (|T^) momentum distributions for the transitions to the 
ground state (circles) and first excited state (squares) in the reaction 7 Li(e, e'p) e He, 
compared to CDWIA calculations with the VMC wave functions (solid lines) for 
AV18/UIX. 



distorted- wave impulse approximation (DWIA) . An analysis of meson factory data 
for p-shell nuclei based on the Cohen-Kurath shell model (|5|) showed that reason- 
able agreement with data could be achieved if the quadrupole excitation operator 
was enhanced by a factor ~ 1.75-2.5 ([[|). Figure [9| shows a calculation of the 
differential cross sections for 7 Li(7r,7r') scattering to the first three excited states 
using VMC transition densities as input, with no enhancement factors (14). Solid 
red lines show the full results, which are in good agreement with the data, while 
the dashed blue lines give the contribution coming from protons only. Because of 
the strong isospin dependence of the ttN scattering i-matrix, this reaction can be 
a stringent test of the neutron densities predicted by the quantum Monte Carlo 
calculations. 

The VMC wave functions for the AV18/UIX model have also been used to cal- 
culate single-nucleon momentum distributions in many nuclei, and a variety of 
cluster-cluster overlap wave functions, such as (dp\t), (dd\ct), and (ad| 6 Li) ([79], 30). 
Recently the overlaps ( 6 He(J 7r ) + p(£j)\ 7 Li) for all possible p-shell states in 6 He 
were studied (13). The spectroscopic factors obtained from these overlaps are 0.41 
to the ground state of 6 He and 0.19 to the 2 + first excited state. These factors are 
significantly smaller than the predictions of the Cohen-Kurath shell model (|8ll), 




Figure 11: Calculated density contours of 8 Be in the lab frame (left) and the 
intrinsic frame (right), labeled with densities in fm~ 3 



which gives values of 0.59 and 0.40, respectively. The CK shell model requires 
that the possible 6 He+p states sum to unity within the p-shell, whereas in the 
VMC calculation, the correlations in the wave function push significant strength 
to higher momenta that cannot be represented as a 6 He state plus p-wave proton. 
The VMC overlaps were used as input to a Coulomb DWIA analysis of recent 
7 Li(e, e'pfRe data taken at NIKHEF @ and found to give an excellent fit to the 
data, as shown in Fig. 1C. 

The VMC wave functions are based on one-body parts that have a shell-model 
structure, namely four nucleons in an a core coupled to (A — 4) one-body £ = 1 
wave functions. However the low-lying states of 8 Be exhibit a rotational spectrum 
and are believed to be well approximated as two a's rotating around their common 
center of mass. It is possible to recover this picture from the VMC wave functions 
by a modified Monte Carlo density calculation (|T9|). 

The standard Monte Carlo method for computing one-body densities, p(v), is to 
make a random walk that samples |\P(ri, r2, • • • , r,4)| 2 and to bin ri, r2, • • • ,ta for 
each configuration in the walk. The density is then proportional to the number of 
samples in each bin. In the case of a J = nucleus, this "laboratory" density will 
necessarily be spherically symmetric. 

The intrinsic density in body-fixed coordinates can be approximated by comput- 
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ing the moment of inertia matrix, A4, of the A positions for each configuration. 
The eigenvalues and eigenvectors of M, are found and a rotation to those prin- 
cipal axes is made. The resulting r^, r' 2 , ■ ■ ■ , r' A is then binned. The eigenvector 
with the largest eigenvalue is chosen as the z' axis. This procedure will not pro- 
duce a spherically symmetric distribution, even if there is no underlying deformed 
structure, because almost every random configuration will have principal axes of 
different lengths and the rotation will always orient the longest principal axis in 
the z' direction. However no artificial structure is introduced flT9|). 

When the above procedure is applied to the three lowest 8 Be states, a dramatic 
intrinsic structure is revealed, as shown in Fig. [ll| for the ground state. The figure 
shows contours of constant density plotted in cylindrical coordinates. The left side 
of the figure shows the standard, lab frame, density calculation. For the J = 
ground state, this is spherically symmetric as shown. The right side of the figure 
shows the intrinsic density. It is clear that the intrinsic density has two peaks, with 
the neck between them having only one-third the peak density; we regard these as 
two a's. This assignment is strengthened by making the same construction for the 
J = 2 + and 4 + states; although the laboratory densities for these states (in M = J 
states) are quite different, the intrinsic densities are, within statistical errors, the 
same as the J = intrinsic density (|i~9|). 

If the + , 2 + , and 4 + states are generated by rotations of a common deformed 
structure, then their electromagnetic moments and transition strengths should all 
be related to the intrinsic moments which can be computed by integrating over 
the projected body-fixed densities. This is explored in Ref. ( |l~9| ) and a generally 
consistent picture is shown. 



7 LOW-ENERGY NUCLEON-NUCLEUS SCATTERING 

The calculations presented so far have treated resonant nuclear states as if they 
were particle stable; that is the VMC trial wave functions decay exponentially 
at large distances and the GFMC propagation does not impose a scattering-state 
boundary condition. As was shown in Fig. |, this GFMC propagation converges 
for unstable states that are narrow, but for wide states, such as 8 Be(4 + ), the 
E(t) keeps decreasing for increasing r. Such states should be computed using 
scattering-wave boundary conditions. 

The 5 He(i", |") states have been studied in VMC (ph and GFMC ph using 



.2 ' 2 

scattering- wave boundary conditions. The VMC calculations use a trial function 
that is a straight-forward generalization of the ^>y given in Sec. 3. In it the ^fj 
contains one p-wave single-particle function that goes to zero at a specified (large) 
radius, R n . The variational energy is minimized subject to this boundary condition. 
The phase shift at the resulting scattering energy, E s = E( 5 He(R n )) — i?( 4 He) is 
then obtained from 

= < 56 > 
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where k = (2/j,E s ) 1 / 2 /%, \i is the a+n reduced mass, and ji and n; are the spherical 
Bessel functions (in this case 1 = 1). It is possible to generalize this formulation 
by specifying a logarithmic-derivative at R n , rather than having the wave function 
go to zero there. This method requires the computation of the energy difference, 
E s . If this is done with two different VMC calculations, one for 5 He with the 
specified boundary condition, and one for 4 He, then the statistical errors of these 
two calculations must be added in quadrature. Instead one can directly evaluate 
the energy difference in a VMC random walk that is controlled by the 5 He wave 
function. In practice this gives E s with a smaller statistical error than either of the 



separate statistical errors (82). The GFMC calculations start with the scattering- 
state trial function and use a modified propagator that preserves the boundary 
condition (|83|). Both methods result in a phase shift and energy corresponding to 
a given boundary condition. By repeating the calculations for different boundary 
conditions, the phase shift as a function of energy can be mapped out. 

The VMC calculations of Ref. (g^) used older two- and three-nucleon potentials 
than have been used for the other calculations reported in this review. They repro- 
duced the qualitative features of the experimental a+n | and | phase shifts, 
but, in particular, gave a spin-orbit splitting about 1 MeV too small. The GFMC 
calculations ( j33"l) used a different, but still somewhat old Hamiltonian and obtained 
values that are closer to experiment. These calculations need to be repeated with 
the modern Hamiltonians described in this review, and for other systems such as 
a+d, 6 He+n, and a+a. 

8 ASTROPHYSICAL ELECTROWEAK REACTIONS 



Most of the key nuclear reactions in primordial nucleosynthesis (84) and solar 
neutrino production (B^) involve only the A < 8 nuclei. Many of these are elec- 
troweak capture reactions that are difficult or impossible to measure in the lab- 
oratory. With the high-precision few-body methods, such as PHH and Faddeev, 
or new effective field theory treatments, the reactions to s-shell final states have 
been calculated with unprecedented accuracy in the last few years. These include 
the 1 H(p, e + z^ e ) 2 H ( [86| ) and 3 He(p, e + f e ) 4 He (|87|) weak capture reactions, and the 
1 H(n, 7) 2 H (38), 2 H(re,7) 3 H and 2 H(p, 7) 3 He (|89[) radiative capture reactions. 

Recently, the VMC method has been applied to some of the radiative captures 
to p-shell final states, including 2 H(a,7) 6 Li, 3 H(a,7) 7 Li, and 3 He(a,7) 7 Be (|i~5| , 
16). As discussed in the previous section, many-nucleon scattering states can be 



constructed within the quantum Monte Carlo framework, but it has not yet been 
done for A > 6 nuclei. Thus for these first studies, the appropriate electromagnetic 
matrix elements (primarily El and E2) are evaluated between a VMC six- or seven- 
nucleon final state, and a correlated cluster-cluster scattering state which is not 
variationally improved. 

Since much of the capture takes place at long range, an important ingredient 
in these calculations is an asymptotically-correct description of the six- or seven- 
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body final state. This feature has to be implemented in the initial variational 
wave function because, as remarked before, the VMC and GFMC methods find 
their best wave functions by optimizing the energy, which is not sensitive to the 
long-range behavior. The asymptotic behavior can be imposed by requiring the 



one-body wave function in Eq.(33) to satisfy the condition: 

[0^[n] (r ^ K Wkm(2ir )/r , (57) 

where Wkm{^l r ) is the Whittaker function for bound-state wave functions in a 
Coulomb potential and n is the number of p-shell nucleons. Here = 2^ n Bi n /fi 2 , 
with fi^n and B^ n the appropriate two-cluster effective mass and binding energy. 
For 6 Li, B42 = 1-47 MeV; for 7 Li or 7 Be, capture can be to the | ground or \ 
excited state with corresponding values for S43. 

The initial-state wave functions are taken as elastic-scattering states of the form 

\i> aT -LSJM) = AU^(r ar )Y LML (r ar )l[G ij \ij^ s ) 1 , (58) 

I lJ ) LSJM 

where curly braces indicate angular momentum coupling, A antisymmetrizes be- 
tween clusters, ip a is the 4 He ground state, and ip™ s is the deuteron or trinucleon 
ground state in spin orientation ms- The Gij are identity operators if the nucleons 
i and j are in the same cluster, else, they are a set of pair correlation operators, in- 
cluding both central and spin-isospin dependent terms, which introduce distortions 
in each cluster, under the influence of individual nucleons from the other cluster. 
They are similar to the correlations discussed in Sec. 3, except that they revert to 
the identity operator at pair separations beyond about 2 fm. The correlations (fi^r 
are generated from optical potentials that describe the experimental phase shifts in 
cluster-cluster scattering; see Refs. (|Tq, 16) for details. Because the Gij go to unity, 



the tp aT has the same phase shifts as those generated by the optical potential. 

Evaluation of electromagnetic matrix elements between the initial scattering and 
final bound states can be split into two parts by noting that all of the energy depen- 
dence is contained in the relative wave function cj)^ and the transition operators. 
Thus, using a technique first applied in a VMC calculation of d+d radiative cap- 
ture (|90|), the matrix element for a given scattering partial wave can be written 
as 



Tf SJiJf (q) = / dxx 2 ^(x) (59) 
Jo 

x <Vd /m/ \Ta(q)A J 8(x - r ar )Y^(r aT )) ]J G^aW) 



LSJiMi 



where Ti\ denotes any of the standard Eg or Mi operators. The integration over 
all coordinates except x can be calculated just once for each partial wave by Monte 
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Figure 12: The total S'-factor for a+d radiative capture calculated for two dif- 
ferent optical potentials with the AV18/UIX nuclear wave functions, compared to 
experimental data (9^,9^,53). 



Carlo sampling, and the result can then be used to compute the full integral for as 
many energies as desired by recomputing <p^(x) only. 

The result for the 2 H(a,7) 6 Li capture 5-factor is shown in Fig. where it is 
compared to a collection of direct and indirect data. The S'-factor falls rapidly 
at low-energy because a pseudo-orthogonality between the initial and final states 
suppresses normal S-wave capture. Whenever there is such a suppression, it is im- 
portant to consider higher-order terms that might contribute, including relativistic 
corrections and two-body charge and current operators (f|). This is possible with 
the VMC calculation because of the fully correlated ^4-body wave functions that 
are used. In this case, a relativistic center-of-energy correction leads to an El con- 
tribution that dominates at low energy, while the resonance region and above is 
primarily E2. Results for 3 H(a,7) 7 Li, and 3 He(a,7) 7 Be can be found in Ref. (16) 

This kind of VMC calculation is only a first investigation into the p-shell astro- 
physical reactions, and should be improved in the future by using the more exact 
GFMC wave functions for both bound and scattering states, and the Hamiltonians 
with improved NNN potentials. 
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9 NEUTRON DROPS 

Neutron-rich nuclei are interesting both because of their importance in various 
astrophysical contexts, such as the r-process or neutron-star crusts, and the cur- 
rent interest in experiments with radioactive beams. These nuclei are often studied 
within the framework of mean-field models using Skyrme or other potential models. 
The parameters of these models are fit to experimentally known binding energies, 
that is for situations with N^Z. In particular the isospin dependence of the spin- 
orbit component of such potentials is considered to be not strongly constrained. 
Neutron drops offer the possibility of theoretical guidance for the isotopic depen- 
dence of such parameters. 

Neutron drops are systems of interacting neutrons confined in an artificial exter- 
nal well. Eight neutrons form a closed shell and single-particle spin-orbit splittings 
can be studied in drops of seven or nine neutrons. Pairing energies can be studied 
if six-neutron drops are also computed. Calculations of systems of seven and eight 
neutrons interacting with AV18/UIX were used as a basis for comparing Skyrme 
models of neutron-rich systems with microscopic calculations based on realistic 
interactions (|39|) . The external one-body well used is a Woods-Saxon: 

v i(r) = T -i r-^r rr^ ; (6°) 

the parameters are Vo = —20 MeV, ro = 3.0 fm, and ao = 0.65 fm. Neither 
the external well nor the total internal potential (vij + Vijk) are individually at- 
tractive enough to produce bound states of seven or eight neutrons; however the 
combination does produce binding. 

Tables [|-|| show results for the neutron drops with this external well. The T = | 
nature of the S^S^. term of V Sn results in large contributions in the neutron drops. 
As a result the seven- neutron drops computed with AV18/IL2 have double the spin- 
orbit splitting predicted by AV18/UIX. Thus the spin-orbit splitting in neutron-rich 
systems depends strongly on the Hamiltonian used. 

10 CONCLUSIONS AND OUTLOOK 

Quantum Monte Carlo methods can now be used to obtain accurate (within 2%) 
energies for light p-shell nuclei up to A = 8. Initial calculations of A=9, 10 nuclei 
have been made by the authors, and studies of ^4=11, 12 nuclei should be feasible 
in the next few years by variational and Green's function Monte Carlo methods. 
This progress in the nuclear many-body problem is due both to the rapid growth 
of computational power and the continuing evolution of algorithms. The new 
auxiliary-field diffusion Monte Carlo method, which samples spins and isospins by 
auxiliary fields, and space by standard diffusion Monte Carlo, may be the key to 
doing even larger nuclear systems (p0[). 
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Studies of the p-shell nuclei allow us to test nuclear forces in new ways not ac- 
cessible in s-shell nuclei, particularly the odd partial waves of NN scattering and 
the T = | triples for NNN forces. They also give us many more cases in which 
to examine charge-dependent and charge-symmetry-breaking interactions. With 
these calculations of nuclear spectra, we see for the first time that nuclear struc- 
ture, including both single-particle and clustering aspects, really can be explained 
directly in terms of bare nuclear forces that fit NN data. A crucial ingredient for 
quantitative agreement is the addition of realistic NNN forces, including at least 
two terms beyond the standard long-range two-pion-exchange potential (22). 

Many aspects of nuclear structure and reactions are described quantitatively in 
these studies. The energy differences between isobaric analog states are explained 
well once the complete electromagnetic and strong charge-independence-breaking 
forces deduced from NN scattering are included. Charge radii and quadrupole 
moments agree with experiment, and we expect magnetic moments will also once 
the important two-body exchange currents are included, as they have been in s-shell 
nuclei (|65|). Elastic and transition form factors measured in electron scattering, 
and transition densities that are tested in pion scattering, are accurately predicted 
without the introduction of effective charges. Spectroscopic factors are naturally 
quenched by the correlations in A-body wave functions and we see that those 
correlations build up intrinsic cluster structure, as in the case of 8 Be. 

Initial calculations of nucleon-nucleus scattering states and electroweak capture 
reactions are encouraging. These two related problems will be a major activity 
for the quantum Monte Carlo studies over the next several years, with studies of 
scattering states in nuclei like 7 He and 9 He, and of the cluster-cluster scattering 
states that enter into the astrophysical capture reactions. Additional reactions like 
7 Be(p, 7) 8 B and even 8 Be(a, 7) 12 C should become feasible. Still other projects will 
include the study of weak decays and, going beyond the p-shell, looking at the 
unnatural-parity intruder states which start to dip into the low-energy spectrum 
with the ^4 = 9 nuclei, and can be particle stable by A = 10. 

Refining the nuclear Hamiltonian will remain a major aspect of future work. 
Extensive new NN scattering data taken since the Nijmegen partial-wave analysis 
has increased the database to more than 6000 points, and a new CD-Bonn 2000 
potential has been constructed that again achieves a x 2 /datum 1 (|94l) . Also, 
intensive searches are being made to look for NNN force signatures in N + d 
scattering experiments (95,96), which may well indicate a need for additional terms 
beyond those contained in the Illinois models. Our ability to test new force models, 
as they become available, in light nuclei by accurate quantum Monte Carlo methods 
will continue to be an important tool for nuclear physics. 
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